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Abstract 



The effects of internal quark structure of baryons on the composition and 
structure of neutron star matter with hyperons are investigated in the quark- 
meson coupling (QMC) model. The QMC model is based on mean-field de- 
scription of nonoverlapping spherical bags bound by self-consistent exchange 
of scalar and vector mesons. The predictions of this model are compared 
with quantum hadrodynamic (QHD) model calibrated to reproduce identical 
nuclear matter saturation properties. By employing a density dependent bag 
constant through direct coupling to the scalar field, the QMC model is found 
to exhibit identical properties as QHD near saturation density. Furthermore, 
this modified QMC model provides well-behaved and continuous solutions at 
high densities relevant to the core of neutron stars. Two additional strange 
mesons are introduced which couple only to the strange quark in the QMC 
model and to the hyperons in the QHD model. The constitution and struc- 
ture of stars with hyperons in the QMC and QHD models reveal interesting 
differences. This suggests the importance of quark structure effects in the 
baryons at high densities. 

PACS number(s): 26.60.+C, 21.65.+f, 12.39.Ba, 24.85,+p 
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I. INTRODUCTION 



Neutron stars are born in the aftermath of supernova explosions with interior tempera- 
tures T ~ 10 11 K, but cool rapidly in a few seconds by deleptonization []J to almost cold 
nuclear matter. Neutron star matter is charge neutral so that gravitational force can bound 
it against the relatively strong repulsive Coulomb force, and is in /3-equilibrium condition, 
i.e. in its lowest energy state. Since the matter density in the core could exceed a few times 
the normal nuclear matter density, neutron star matter provides an interesting possibility 
to investigate the strong interaction effects which are poorly understood at supernuclear 
density. In fact, the structure of a neutron star is chiefly determined by the equation of 
state (EOS) of the strongly interacting constituents. 

There have been several attempts to determine the EOS for dense nucleon matter which 
are based primarily on nonrelativistic potential models and relativistic field theoretical mod- 
els. The former approach comprises a Hamiltonian with a two-nucleon potential fitted to 
nucleon-nucleon scattering data and the properties of deuteron. The quantum many-body 
problem is traditionally handled either by a selective summation of diagrams in perturba- 
tion theory (the Brueckner-Bethe-Goldstone approach) or by using a variational method 
with correlation operators ||||. However, an important shortcoming of many potential 
models are that they are well suited at low densities only; the EOS becomes acausal i.e., 
the speed of sound exceeds that of light at densities relevant for maximum mass neutron 
stars. Moreover, these models lead to symmetry energies that drastically decreases beyond 
about three times the saturation density which is a serious deficiency for highly asymmetric 
systems like neutron stars. The alternative approach to the nuclear many body problem 
involves the formulation of effective relativistic field theory within the framework of quan- 
tum hadrodynamics (QHD) where the appropriate degrees of freedom are the baryons 
interacting through the exchange of isoscalar scalar and vector mesons (a, uj) and the isovec- 
tor vector p meson which provides the driving force to the isospin symmetry. The theory 
is an effective one since the coupling constants are determined by the saturation proper- 
ties of nuclear matter. The equations of motion for the baryons and mesons are solved 
self-consistently in the mean field approximation. There exists a large body of calculations 
based on QHD which were found to provide a realistic description of bulk properties of finite 
nuclei and nuclear matter 0-0. A large number of single-particle properties of finite nuclei 
are also well accounted including the charge distributions, spin-orbit interactions etc 



The central to the success of the theory is that the small binding energy in a nucleus arises 
from the cancelation between large Lorentz scalar and vector potentials, each of which is 
approximately several hundred MeV at saturation density. The potentials being compara- 
ble to the nucleon mass, it indicates the importance of relativistic effects even at normal 
densities. With increasing density the Fermi momentum of the baryons increases further 
which points to stronger relativistic effects. Therefore QHD appears to be more appropriate 
than nonrelativistic models for neutron star calculations. Moreover, the Lorentz covariance 
in this theory is retained from the outset. Consequently, the EOS automatically respect 
the causality limit. The relativistic mean field models with baryons and mesons as point 
particles, and the coupling constants determined from the saturation properties of nuclear 
matter, has been extrapolated into the high density regime to investigate the neutron star 



properties 11-15 
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At high densities in the core of a neutron star the relevant degrees of freedom that 
may be crucial to its structure and composition are the effect of quarks confined within 
the hadrons. Small but interesting corrections to the standard hadronic picture has been 
already observed such as the EMC effect which reveals the medium modification of the 
internal structure of the nucleon ||16|| . By treating the nucleons as structureless particles, 
the QHD model completely misses the important effect of the constituent quarks especially at 
large densities. It is therefore instructive to study the relevance of the quark structure of the 
nucleon at supernuclear density regime as in the neutron stars. Quantum chromodynamics 
(QCD) which governs the underlying strong interactions of quarks and gluons, although, 
believed to be essentially the fundamental theory at the nuclear and subnuclear scale is 
however intractable due to the nonperturbative features of QCD. It is therefore not a theory 
from which one can derive practical results for the equation of state. The most fruitful way 
to deal this situation would be to approach the problem from both sides, i.e. to develop the 
QHD and QCD motivated models and to compare their predictions. 

The spontaneous symmetry breaking and chiral symmetry restoration criteria have been 
utilized to model effective Lagrangians for the low energy strong interaction, as for exam- 



ple the Nambu— Jona-Lasinio model [[L7|]. However, these models have never been applied 



successfully for the description of the saturation properties of nuclear matter. On the other 
hand, the quark-meson coupling (QMC) model proposed by Guichon [18] provides a simple 
and attractive framework to investigate the direct quark effects in nuclei. The model de- 
scribes nuclear matter as nonoverlapping, spherical and static MIT bags in which the quarks 
interact through a self-consistent exchange of structureless scalar a and vector uo mesons in 
the mean field approximation. This simple QMC model has been subsequently refined by 



including nucleon Fermi motion and center of mass corrections to the bag energy |T9[ and 
applied with reasonable success to various problems of nuclear matter ppj-pB]] and finite 
nuclei [f26| -|28"11 ■ Recently, the model has also been used to investigate the properties of A, £ 
and H hypernuclei |f2U| . 

Although it provides a simple and interesting framework to describe the basic features 
of the nuclear systems in term of quark degrees freedom, the QMC model has two serious 
shortcomings. Firstly, it predicts much smaller scalar and vector potentials for the nucleon 
than that obtained in the well-established QHD model. This implies a much smaller nu- 
cleon spin-orbit potential (sum of scalar and vector potentials) and hence weaker spin-orbit 
splittings in finite nuclei. Secondly, it is not apparent how far the EOS in the QMC model 
can be extrapolated at high densities in the neutron star interior because the assumption 
of nonoverlapping bags may break down. It was recently pointed out |24|,|25|] that the small 
values of the nucleon potentials stem from the assumption that the bag constant is fixed 
at its free space value. By introducing a density dependent bag constant (by considering it 
to be proportional to higher powers of the scalar field a, for example) it was demonstrated 



25| that large scalar and vector potentials are obtained. Evidently, a drop in the bag value 



with increasing density implies a decrease in the effective nucleon mass in a self-consistent 
manner. Furthermore, it was shown that the quark substructure is entirely contained in the 
scalar potential |20[| . Therefore, the density dependent bag constant through a field repre- 
sents, in a way, self-interaction in the scalar field in this model. This feature of nonlinear 
interaction in sigma field, obtained by a field dependent bag constant may also remedy the 
other shortcoming of the overlapping bags at high density matter. One possible way to verify 
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that the QMC model may be applied at supernuclear densities is by exploring observables as 
a function of density which shall not deviate from its continuous behavior at high densities. 
In other words, the QMC theory which provides a reasonable behavior of the equation of 
state etc. in the vicinity of saturation, should also be well-behaved when extrapolated to 
supernuclear density regime. Of course, the QMC model at densities near nuclear matter 
value should reproduce the results based on the more established QHD model. Any devi- 
ation (not discontinuity) in the QMC results from that of the QHD at high densities may 
then be regarded as the effects of the quark structure beyond the structureless hadronic 
picture. 

In this paper we shall show that by employing a bag constant coupled to the a field, the 
continuity of the equation of state etc. indeed persists in the QMC model when extrapolated 
to high densities relevant to the neutron star core. This will allow us to investigate the 
structure and composition of a neutron star in the QMC model and compare its predictions 
with that in the QHD model. Since the QMC model with bag constant as function of a 
field represents nonlinear scalar self-interactions, for comparison we employ a QHD model 
which contains cubic and quartic scalar self- interactions PD| ; the latter model is calibrated 
to reproduce the same saturation properties as the QMC model. 

At high densities and hence large nucleon Fermi energies expected in the cores of stars, 
weak interaction energetically favor the conversion of some nucleons at the Fermi surface 
to hyperons. The negatively charged hyperons can then also replace the leptons. These 
effects cumulatively can decrease the pressure resulting in the softening of the equation 
of state. Within the QMC framework, only hypernuclear matter and hypernuclei |]22|j29[ 



(without charge neutrality and /^-equilibrium conditions) have been investigated where the 



bag constant is fixed to its free space value. Moreover, it was assumed p2| , |29|l that the 
nonstrange a and u> mesons couple only to the u and d quarks; the s quark is unaffected in 
the medium and set to its constant bare mass value. We investigate neutron star properties 
with an improved Lagrangian by incorporating an additional pair of hidden strange meson 
fields (a*, <ft) which couple only to the s quark in the QMC model and only to the hyperons 
in the QHD model. The standard mean field model was found to be inadequate to describe 



the strongly attractive hyperon-hyperon interaction observed in double A hypernuclei [31 



The hyperon-hyperon interaction is expected to be important for hyperon rich matter flT 
present in the cores of neutron stars. 

The paper is organized as follows. In section II, we introduce the relativistic mean 
field models, the QMC and the QHD, both extended to incorporate the two additional 
(hidden) strange mesons. The relevant equations for neutron star matter with hyperons are 
summarized in these models. In section III, we make a systematic comparison of the QMC 
and the QHD results in the entire density range relevant to neutron stars. In section IV, 
we use the model to study the constitution and structure of a star. Section V is devoted to 
summary and conclusions. 



II. THE RELATIVISTIC MEAN FIELD MODELS 
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A. The Extended Quark-Meson Coupling Model 



In this section we present a brief introduction to the quark meson coupling model for 
baryonic matter and its extension to include the hidden strange mesons which couple explic- 
itly only to the s quark in a hyperon bag. In the QMC model, a baryon in nuclear medium 
is assumed to be a static spherical MIT bag in which the quarks are coupled to the me- 
son fields in the relativistic mean field (RMF) approximation. In QMC model calculations 
for hadronic matter 



22,29 



the meson fields considered are isoscalar scalar a and vector u 
mesons and the isovector vector p 5 " meson. However, they being nonstrange are allowed to 
couple only to the u and d quarks within a baryon bag, while the s quark is unaffected. It 
is expected that with increasing density, the s quark mass should also be modified. To im- 
plement this situation, we incorporate two additional mesons, the scalar meson fo(m = 975 
MeV) (denoted as a* hereafter) and the vector meson <p(m = 1020 MeV) with their masses 
given in the parenthesis. These strange mesons couple only to the s quark in a hyperon bag. 
This extended QMC model has an additional advantage that it accounts for the strongly 
attractive AA interaction observed in hypernuclei which cannot be reproduced by (a, lu, p) 
mesons only |3T . 



For a uniform static matter within the RMF, let the mean fields be denoted by a, a* for 
the scalar mesons, and u Q , O an d P03 f° r the timelike and the isospin 3-component of the 
vector and the vector-isovector mesons. The Dirac equation for a quark field ip q in a bag is 
then given by 



27 • d 



9> 



q * 
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9luo + 9l<Po + 2#p r ^Po3 



(1) 



where g%, g**, g« 



g q are the quark-meson coupling constants and m q the bare mass of 



the quark q = (u, d, s). The normalized ground state for a quark in a bag is given by 



ip q (r, t) = Af q exp (-ie q t/R 



Bj 



jo (x q r/R B ) 
i/3 g a ■ rji (x q r/R B ) 



Xg 



(2) 



where 



Q q + Rb( 9l^o + 9l<Po + 2#p r *Po3 



0o 



fi g - R B m* q 
^tt q + R B m* 



(3) 



The normalization factor is given by 



K 2 



2R B j 2 (x q ) n q (n q -i) + R Bm ;/2 



x„ 



(4) 



with Q q = x q + (R B m*) 2 the kinetic energy of the quark q and R B is the radius of a baryon 
B, and \ q the quark spinor. The effective mass of a quark is given by 



m q - g H a a 



q * 



(5) 



The linear boundary condition, jo(x q ) = (3 q ji(x q ), at the bag surface determines the eigen- 



value x q . 
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The hyperon couplings are not relevant to the ground state properties of nuclear matter, 
but information about them can be available from levels in A-hypernuclei ||32|| . Experimental 
data of E-hypernuclei are scarce and ambiguous because of the strong EiV — > AN decay, 
while only few events in emulsion experiments with K~ beams have been attributed to the 
formation of H~ hypernuclei. Therefore considerable uncertainty in the hyperon-nucleon 
interaction exists even at the normal nuclear density, their interactions at high densities are 
more ambiguous. In view of the uncertainties in the hyperon couplings, for simplicity, we 
employ in this paper the SU(6) symmetry based on the light (u, d) quark counting rule for 
both the scalar (a, a*) and vector (uj, <fi) coupling constants to the hyperons. The coupling 
constants are thus related by 

111 

^9aN = -jQak = -j9u-il = 9aE = 9a , 
~9wN — ~9ljA — ~9ljT, — QuiE = 9u; , 

9 p n = ^9 p t, = g P s = 9p , 9 P k = . (6) 
The couplings to the strange mesons are 

2\p2 r 
2g a * A = 2g a *x = g a * E = -^~9aN = 2V2g q cr , g a * N = 

2^/2 

2^a = 2^ = g^ = -^-g^N = 2V2gl , g^ N = . (7) 

Note that in a baryon the u and d quarks are not coupled to the strange (a*, 0) mesons 
i.e., g^'* d = g^ ,d = 0, while the s quark is unaffected by the (a, u) mesons i.e., g s a = g^ = 0. 
Furthermore, from Eqs. (6) and (7) it is evident that the coupling constants of the s quark 
to (<j*, 0) may be obtained from the coupling constants of the (u, d) quarks to (a, uj) mesons 
by the relation g s a * = v2g%' d and g^ = \[2g^ d . The energy of a baryon bag consisting of 
three ground state quarks is then given by 

E%* = + *& bBb , (8) 

Kb o 

where n q is the number of quarks of type q, zb accounts for the zero-point motion, and Bb 
is the bag constant for the baryon species B. After the corrections of spurious center of 
mass motion, the effective mass of a baryon is given by [1£,2C 



m B 



(< g ) 2 -E^(V^) 2 . (9) 



For fixed meson fields, the bag radius Rb is determined by the equilibrium condition of the 
baryon bag in the medium dm* B / 'ORb = 0. For a given value of bag constant Bb = B in 
free space, the parameter zb and the bag radius Rb of the baryons may be obtained by 
reproducing the physical masses of the baryons i.e. Eq. (9) in free space. In the present 
calculation, the current quark masses considered are m u = rrid = and m s = 150 MeV. 
For our choice of free space bag constant B^ 4 = 188.1 MeV, the values of zb and Rb are 
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collected in Table I. It is seen that for the fixed bag value, the equilibrium condition in free 
space results in an increase of the bag radius and a decrease of the zero-point motion for 
the heavier species. 



In the original version of the QMC model [18[ the bag constant was held fixed at its 
free space value B = Bq. The bag constant is a nonuniversal quantity associated with the 
QCD trace anomaly. When a baryonic bag is immersed in matter, it is expected to decrease 
from its free space value as argued in Ref. J53 1 . At present, however, no reliable information 
on the medium dependence of Bb is available on the level of QCD calculations. Effective 
models, as for example the Nambu— Jona-Lasinio (NJL) model |T7[ which approximate low 
energy QCD are constructed based on symmetries and symmetry breaking patterns of QCD, 
in particular, the chiral symmetry breaking. The concept of bag constant arises naturally 
in these models where its value decreases when the density of the nuclear environment is 
increased [34|. To reflect this physics in the QMC model, the density dependence of the 
bag constant for a nucleon was proposed by Jin and Jennings by considering its direct 
coupling to the scalar mean field, i.e., 



B N /B = exp -Ag^a/ m N 



(10) 



where is a real parameter; in this paper we shall adopt this form of exponential depen- 
dence. This direct coupling model is inspired by NJL type nontopological soliton model for 



the nucleon |35| where a scalar soliton field is responsible for the confinement of the quarks 
to form a nucleon. When the nucleon soliton is inserted into the nuclear environment, the 
scalar soliton field will interact with the scalar mean field [36|. In fact, a similar approach 



is used to construct in QHD the nonlinear mean field models, where the unknown density 
dependence of the nuclear energy functional is parametrized by nonlinear meson-meson in- 
teractions |37|j38[1 . The direct coupling of the bag constant to the scalar mean field a in 
nucleonic medium needs to be extended for hyperons where additional scalar field a* is em- 
ployed. To this end, we employ consistently the SU(6) symmetries for the scalar couplings 
of Eqs. (6) and (7). Specifically, the bag constant Bb of a baryon is directly coupled to o 
and a* fields through the relation 



B B /B = exp 



{ £ n q a + (3 

\q=u,d 




m B 



(11) 



where m# is the bare mass of the baryon B, and for nucleonic bags n q — 3. This modeling 
has only one real positive parameter g® which for nucleon is related to of Eq. (10) by 
gj* = g^ /3. We have refrained from using an additional parameter for coupling to the a* and 
have rather used the SU(6) symmetries because we believe that a reliable extrapolation to 
high densities should be based on a model having as few adjustable parameters as possible so 
that the model having been fitted to saturation properties of nuclear matter can be tested for 
its predictive power under conditions not included in the determination of the parameters. 
It may be also noted that in the parametrization (11), the use of free baryonic mass mg 
is essential. By considering only (er, ou, p) mesons in the model, we have found that this 
choice of B B /B results at nuclear matter density the scaling relation 5m* AS /5m* N « 2/3 
and 5mZ/Sm* N ~ 1/3, where 5m* B = rriB — rn* B for the baryon B. If we now define the field 
dependent a — B coupling constant ^ CT s(cr) by 
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g a B(cr)cr = m B - m* B (a) , 



(12) 



the same scaling relation is obtained at the nuclear matter density i.e., 5m* B /Sm* N = 
9<jb{o) / 9aN{o) are 2/3 and 1/3 for (A, S) and S, respectively. This implies that the 
parametrization (11) in the nonstrange sector is consistent with the SU(6) symmetry em- 
ployed in determining the a — B coupling constants. In this model the baryon effective mass 
m* B and the bag constant are determined self-consistently by combining Eqs. (8), (9) and 
(11). In principle, the parameter zb may also be modified in baryonic medium. However, 
unlike the bag constant it is not clear how zb changes with density as it is not directly 
related to chiral symmetry. In this paper we assume zb remains constant at its free space 
value zb = z as given in Table I. 

The total Lagrangian density of a neutron star matter for the full baryon octet in the 
QMC model within RMF approximation can be written as 

£qmc = E h ■ d - m* B (a,a*) - 7 {g^s^o + g^o + t^b^Pc^ 4>b 

+ \ ( m W + ml*a* 2 + mlul + m\^\ + m 2 p p 2 03 ) + (17 • - m t ) ip t . (13) 

Here the sum on B is over all the charge states of the baryon octet (p, n, A, S + , E°, £~ , S°, S _ ) 
coupled to the a,u, p,a*,(f) mesons. The sum on I is over the free electrons and muons (e - 
and p~) in the star. The baryon effective mass of Eq. (9) may be expressed in terms of the 
field dependent a — B and a* — B coupling strengths g a B(&) and g a *B(&*) as 

m* B (a, a*) = m B - # ct b(<t)<t - g a * B (cr*)a* . (14) 

The dependences of the coupling strengths on the applied scalar fields must be calculated 
self-consistently within the quark model. 

For the QMC model, the equations of motion for the meson fields in uniform static 
matter are given by 



m 2 c a 



^ . ,2J B + 1 f k B m* R (a,a*) l2 ,, . . 



2 _* ^„ r < 1 * ^ J b + 1 f kB rn* B (a,a*) 2 

" g-Btoi* J q ——^ kdkt (16) 

™>o = £<7u,b (2Jb + 1) 4/(6vr 2 ) , (17) 



B 



m p 



mlfo = E 9<t>B (2J B + 1) 4/(6tt 2 ) , (18) 

B 

P03 = E 9 P BhB (2J B + 1) k 3 B /(6n 2 ) . (19) 
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In the above equations Jb and I%b are the spin and the isospin projection and ks is the 
Fermi momentum of the baryon species B. On the right hand sides of Eqs. (15) and (16), 
a new characteristic feature of QMC beyond QHD appears through the factors C B {cr) and 
Cb(c*) where 



9aBC B {cr) 




1 - 



fin 



Sb(<?) + 



E b B as 



m* B (a,a*) 3 



(20) 



g** B C B (a ) = w~ = n s gl* —— , 1 - . - S B (a ) + 



da* S3a m* B (a,a*) \\ E h B ^R B ) V ' E h B ^ j 

+ , (21) 

m B (<j, a*) 3 m B 

where n s (= 3 — J2q= u ,d n g) * s ^ ne number of s quark in the baryon. The quark scalar densities 
in the bag are 

, x r - n q /2 + R B m*(n q -i) 

w°1*^ = 4to^ ! (22) 

and a similar expression for Sb (&*) for the contribution from the medium modification of 
s quark in the field a*. The medium dependence of the scalar densities on the bag radius 
was found to be rather insensitive. The last terms in Eqs. (20) and (21) originate from the 
density dependent bag constants through direct coupling to the scalar fields. It is evident 
from Eqs. (11), (15) and (16) that an exponential dependence of the bag constant on the 
fields a and a* introduce a nonlinear self-interaction in these fields. Moreover, the decrease 
of Cb with increasing density provides a new source of attraction and thereby constitutes a 
new saturation mechanism which is different from QHD. 

For stars in which the strongly interacting particles are baryons, the composition is 
determined by the requirements of charge neutrality and /3-equilibrium conditions under the 
weak processes B x — > B 2 + I + V h and B 2 + I — > Bi + v x . Under the conditions that the 
neutrinos have left the system, the charge neutrality condition gives 

(Ztat = E Vb(2Jb + l)4/(6vr 2 ) + £ Ql kf /(3tt 2 ) = , (23) 

B l=e,fi 

where q { corresponds to the electric charge of species i. Since the time scale of a star is 
effectively infinite compared to the weak interaction time scale, weak interaction violate 
strangeness conservation. The strangeness quantum number is therefore not conserved in 
a star and the net strangeness is determined by the condition of /3-equilibrium which for 
baryon B is then given by Hb = b B fi n — qB^ei where \x B is the chemical potential of baryon 
B and b B its baryon number. Thus the chemical potential of any baryon can be obtained 
from the two independent chemical potentials /i n and fi e of neutron and electron. The Fermi 
momentum of the baryons can be obtained from the solution of the equation e B (k B ) = (J>b, 
where the energy eigenvalues of the Dirac equation for the baryons are 
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e B {k) = \Jk 2 +m*g{a,a*) + g uB Wo + g^B&o + 9 P BhBpo3 



(24) 



The lepton Fermi momenta are the positive real solutions of (k 2 + m 2 e ) 1 / 2 = « e and (k 2 + 
m ^) 1 ^ 2 = AV = ^e- The equilibrium composition of the star is obtained by solving the set 
of Eqs. (15)-(19) in conjunction with the charge neutrality condition (23) at a given total 
baryonic density n B = J2b^b(^Jb + 1)^b/(6tt 2 ); the baryon effective masses are obtained 
self-consistently in the bag model. The total energy density and pressure including the 
leptons can be obtained from the grand canonical potential to be 



+ 



1 22,1 2 *2,1 22,1 2 2 

1/2 



2 2 
-, m P P03 



E 



2J B + 1 
2tt 2 



k 2 dk 



E 



JO 



k 2 + mf 



1/2 



/c 2 dk 



(25) 



n ^ 2 2 1 _2 *2 , 1 2 2 i 1 2 j.2 , 1 2 2 

p = "2 m,T(T ~ 2 ^ 2 m< ^° 2 m *^° 2 mpPo3 

1 2 J B + 1 A; 4 dk 1 V 1 /"*' ^ 4 

+ 2tt 2 io [P + mf(a )( T*)] 1/2 + 3Y^^ + " 

To obtain the coupling constants and the parameters in the QMC model, we recall that 
for a fixed B^ 4 = 188.1 MeV, the z values have been adjusted to reproduce the baryon 
masses in free space and these are listed in Table I. For a given value of g%, once the 
three coupling constants g W N, gpN and gj 3 are adjusted, the other coupling constants of 
the hyperons to the meson fields can be obtained by employing the SU(6) symmetry from 
Eqs. (6) and (7). For this purpose, the QMC model is solved for symmetric nuclear matter, 
and as in Ref. [p5|| , for a given value of g q a = 1 the coupling constants g^N, gpN and g^ 



are adjusted to reproduce the nuclear matter binding energy B/A = 16 MeV at saturation 
density no = 0.17 fm -3 and symmetry energy a sym = 32.5 MeV. The resulting coupling 
constants are given in Table II. For the parametrization employed here, the predicted values 
of effective nucleon mass and compressibility at saturation density are m* N /m^ = 0.78 and 
K = 289 MeV. It may be worth mentioning that in order to reproduce the same saturation 
properties of density, binding and symmetry energies, the coupling constants required in 
the original QMC model with bag fixed at E 1/4 = B^ = 188.1 MeV are # 2 /4vr = 20.2, 
g 2 /4:7r = 1.55 and g 2 /^ = 5.51 with a relatively larger effective mass m* N /rriN = 0.89 
and smaller compressibility K = 220 MeV. The higher effective mass and thereby smaller 
scalar field potential is compensated at the saturation density by a smaller vector field i.e., 
a smaller coupling g 2 /ATT. Since at high densities the vector field dominates over the scalar, 
the smaller vector coupling leads to a softer EOS. The parameters obtained here are entirely 
from free space value and from nuclear matter at the saturation density, therefore this set 
can be used also in the model with the two strange mesons. In the following we refer to the 
model where the interaction is mediated by (a, u, p) mesons as QMCI while its extension 
by incorporating (a*, 0) mesons as QMCII. 
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B. The Extended Quantum Hadrodynamics Model 



Quark meson coupling models are designed to describe both the bulk properties of nuclear 
systems and medium modifications of the internal structure of the baryon. Before any 
reliable predictions for changes due to the quark substructure can be made especially at large 
densities relevant to the core of neutron stars, it is important that the QMC model predicts 
the established results of nuclear phenomenology obtained in the quantum hadrodynamics 
model (QHD) near the saturation density. In QHD the relevant degrees of freedom are 
the structureless baryons interacting by the exchange of mesons. The direct coupling of 
the bag constant to the scalar field in the QMC model mimics scalar self-interaction terms. 
Therefore, for a consistent comparison with this QMC model, we employ a version of QHD 
model which contains the cubic and quartic scalar self-interactions 1 30f| - The Lagrangian for 



the baryon octet in the QHD model within the RMF approximation is similar to that of Eq. 
(13) for the QMC model, but for the baryonic effective mass given by 



m 



B 



m B - gaBcr - g a *Bcr* , (27) 



and contains a scalar self-interaction term 



U(a) = fa 3 + fa 4 , (28) 



proposed by Boguta and Bodmer [[!(]] to get a correct compressibility at saturation density. 
Comparing Eqs. (14) and (27), it is clear that the coupling constants in QHD are indepen- 
dent of the scalar field and they are determined at the saturation density. In fact it has 
been demonstrated that QMC model is formally equivalent to the nuclear QHD model 
with a field dependent scalar-nucleon coupling. The equations of motion for only the scalar 
meson fields of Eqs. (15) and (16) are then modified in the QHD model to 

2 d ^ 2J B + 1 r k B m * 



2J B + 1 r k B m 



! - -I>-b^F / „,, - ' „« *» «* • (30) 



B 



2tt 2 Jo [F + mf] 1/2 



Here nonlinear interaction only for the a meson is employed, the interaction between the 
hyperons in the QHD model are through linear a* and 4> mesons. The fact that all the other 
meson field equations are unaltered in QHD, suggests that in the QMC model the effect 
of the internal quark structure of a baryon enters entirely through the factor C B (a) and 
C B (a*) as has been mentioned in Ref. p0| . The vector fields in QMC cause only a shift in 
the quark wave functions. The equations of motions for the meson fields in QHD are solved 
self-consistently in accordance with the charge neutrality and /^-equilibrium conditions to 
obtain the composition and structure of a neutron star. The five coupling constants g a N, fiWv, 
gpN, #2 and #3 in this model are determined by reproducing the same equilibrium properties 
of saturation density, binding energy, symmetry energy, effective mass and compression 
modulus of the QMC model; these are given in Table II. All the other coupling constants 
can be obtained from these couplings. The constant g% in the scalar field potential (28) is 
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found positive. This avoids the fatal problem for a quantum field theory that the energy 
functional may be unbounded from below which leads to instabilities at high densities with 
large scalar fields. Hereafter we refer to the QHD model with (a, u, p) mesons as QHDI 
and its extension by introducing (a*, 0) mesons as QHDII. 



III. RESULTS AND DISCUSSIONS 

In this section we shall present results for baryonic matter in charge neutral and /3- 
equilibrium conditions appropriate for a neutron star in the QMC and QHD models. The 
effective baryon masses m* B /mB defined in Eq. (9) are shown in Fig. 1 as a function of 
baryon density % for the models QMCI and QMCII. Unless otherwise mentioned the thin 
lines refer to results for different species in the model QMCI, while the thick ones correspond 
to those in the QMCII model. The effective masses of the nucleons rapidly decrease with 
increasing density and then saturate at higher densities. Since the nucleons do not couple 
to the strange scalar field a*, their masses in the models I and II are identical. At densities 
around nuclear matter values where hyperons have not appeared, the effective mass of a test 
hyperon is determined by only the scalar field a created by nucleons which is assumed to be 
unaffected by inserting the hyperon. Consequently, the effective masses of all the baryons for 
KB < 2n reveals no splitting in models I and II. Moreover, the different baryonic masses at 
the saturation density indicates the SU(6) symmetry (based on the number of light quarks 
counting rule) for a — B coupling i.e., 5m* AJ] /5m* N = 2/3 and Smi/Sm* N = 1/3. The QMCI 
model respects this scaling relation to nearly the entire density range explored here. On the 
other hand, in the QMCII model at densities n B ~ 0.38 fm~ 3 when the hyperons A) 
production threshold is reached (see Fig. 6), the attractive scalar field a* starts to contribute. 
The reduction in the mass of strange quark in model II entails a substantial decrease in the 
masses for the hyperons in accordance with the SU(6) relation for a* — B coupling. Since the 
S hyperon has two strange quarks, its mass is maximally suppressed by the a* field. Note 
that although m* N /mN reaches small values at high densities in presence of hyperons, it never 
becomes negative in the density range studied here. In fact, relativistic mean field models 
fitted to the bulk properties of nuclear matter with a high m* N /m N ^0.7 — 0.82 lead to a 
finite effective mass even at central densities of maximum mass stars. However, the effective 
masses in this case become negative only at densities much higher than the central densities 



of maximum mass stars [12]. In principle parameter sets fitted to finite nuclei properties 
can also be obtained. However, in most of these sets the effective masses of nucleons in 
presence of hyperons get negative at densities much smaller than that of maximum mass 
stars. Therefore, these sets are not reliable to calculate neutron star properties. Recall that 
the original QMC model with bag constant fixed to the free space value, leads to much 
smaller scalar field and higher effective nucleon mass m^/mjv = 0.89. To reproduce the 
correct spin-orbit splitting, the reduction in the effective mass (m* N /m N = 0.78) is achieved 
by direct coupling of the bag constant to the scalar fields. 

The variation of bag constant for the baryons Bb/Bq (see Eq. (11)) with density is shown 
in Fig. 2 for the QMCI and QMCII models. The bag constants decrease with density and 
saturate at high densities. This behavior is similar to that of effective baryon masses as they 
have been obtained in a self-consistent manner. The decrease of bag constant relative to its 
free space value implies a decrease of bag pressure which causes an increase of bag radius in 
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the medium. The variation of the bag radius Rb/Ro (relative to its free space value Rq which 
is different for different baryons; see Table I) with density is shown in Fig. 3. At saturation 
density when the bag constant for nucleon has decreased to B^/Bq = 0.45, the corresponding 
radius has increased to R^/Ro = 1.22. At densities of n B = (6 — 8)n corresponding to the 
maximum mass of neutron stars, the nucleon bag constant has decreased to a significantly 
small value Bn/B = (0.093 — 0.065) while the corresponding radius is 75 — 88% larger than 
its free space value. This implies a considerably swollen nucleon (and hyperons) in the star 
matter; a detailed discussion of its consequences is given later. 

We now present a systematic comparison between the QMC and the QHD models for 
neutron star matter with their coupling constants determined to reproduce the same set of 
nuclear matter saturation properties as given in Table II. In Fig. 4, the baryon effective 
masses m* B /mB as a function of density n B are displayed for the models QHDII (thin lines) 
and QMCII (thick lines). It is seen that at low and moderate densities, the m* B s for the two 
models are in good agreement, this is not surprising as both the models are calibrated to the 
same properties at nuclear matter densities. At higher densities, especially when hyperons 
start to populate (at n B ~ 0.38 fm~ 3 ), the effective masses are rather distinct in the two 
models. The pure scalar and vector field strengths are shown in Fig. 5 as a function of density 
for the models QMCII (top panel) and QHDII (bottom panel). The potentials for a given 
baryon species are obtained by multiplying them with the corresponding coupling constants 
listed in Table II; for the scalar fields these couplings are however density dependent for 
QMCII. A careful examination of Fig. 5 indicates that the values of the fields and potentials 
for uj in the two models are nearly identical over the entire density range. On the other hand, 
at all densities the a field in QMCII is larger than that in QHDII. However, the decreasing 
coupling constant g a n{a) with n B in the former model causes the potential U a = g a N(cr)a 
to be the same as QHDII at the normal nuclear matter density. This lead to the same 
saturation properties (binding energy and density) in the two models. At densities higher 
than the normal nuclear matter value, go-jv(cr) further decreases causing the potential U a in 
QMCII to saturate earlier than QHDII. In other words, the scalar density factor Cn(<j) [see 
Eqs. (15) and (20)] in QMCII decreases with increasing g q a a (or n B ) as quarks in the medium 
become more relativistic p0| , p9| . As a consequence, the drop in the nucleon effective mass 
relative to its free space value in QMCII is smaller than that in QHDII. Clearly at high 
densities the quark substructure of the nucleon plays a crucial role in QMC model. This 
feature of larger effective masses in QMCII is more evident for the hyperons. This is because 
of a smaller value of the attractive scalar field a* (see Fig. 5) and a decreasing in-medium 
coupling constant g a * B (o-*) results in a much smaller potential U a * = ga*B( cr *) cr * f° r QMCII 
compared to that in QHDII. The effect is most pronounced for S having two strange quarks. 
We have observed similar qualitative differences in models I, i.e. between QMCI and QHDI. 
However, the distinction in m* B for hyperons in the two models is not so profound as they 
lack the a* meson. This indicates the importance of the strange mesons (a*, <fi) which helps 
in revealing more clearly the quark structure of the baryons at high densities in the QMCII 
model as compared to the structureless baryons inherent in the QHDII model. 

Having investigated the crucial role played by quarks confined within the baryons at high 
densities, we shall now consider whether the QMC model based on nonoverlapping bags can 
at all be extended to densities appropriate to neutron star interior. As is evident from 
Fig. 3, at rather large densities the increasing bag radius implies considerable overlapping 
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between the bags and the nonoverlapping bag picture of nuclear matter may break down 
because the effects of short-range correlations among the quarks which should be associated 
with the overlap of hadrons are explicitly neglected in the quark-meson coupling model. 
We however observe that in the present QMC model the physical observables for instance 
the effective masses, fields and therefore the equations of state (see Fig. 7) indicate a 
smooth and continuous behavior from low to very high densities without any dramatic 
discontinuity. Any deviation at the high density regime from the QHD results may be 
interpreted as the interesting effects due to quark structure of the baryons only. This is in 
consonance with the argument put forward j|0| that the physical observables do not depend 
on the bag radius. On contrary, in the original QMC model with bag fixed at the free 



space value, the observables were found |39j to differ drastically from QHD results with 



large discontinuities. Moreover, solutions in that version of QMC model cease to exist above 
ris ~ 4.92n since the eigenvalue x q in Eq. (2) vanishes. This indicates that by using 
a density dependent bag constant through direct coupling to the scalar field, the QMC 
model not only reproduces the correct spin-orbit potential but possibly also includes the 
effects of quark-quark correlations associated with overlapping bags which was missing in the 
original QMC model. A possible explanation to this may follow from an alternative approach 
inspired by effective field theories. In this theory the constraint of renormalizability of the 
hadronic Lagrangian is abandoned |37|j38[| which allows to introduce higher order meson 



self-interaction (i.e. orders higher than the quartic scalar interaction) consistent with the 
underlying symmetries of QCD. By a suitable truncation at some low orders of the fields, it 
was argued [0 that by allowing non-linearities in the meson fields that generate additional 
density dependence in the interactions, the important effects of correlations (and exchange) 
between the nucleons beyond the simple Hartree contributions are automatically included. 
Thus, one can include many body effects beyond the simple Hartree level even though only 
classical meson fields and local interactions are retained. 

In the context of QMC model, as demonstrated in this paper and also in Ref. [BU 



an exponential dependence of bag constant Bb on a (and a*) is equivalent to scalar self- 
interactions of infinite orders, of course each term is smaller by a factor ~ (g^ /tub) than 
its preceding one [see Eq. (11)]. Hence these higher order self-interaction terms possibly 
include (implicitly) the effects of quark-quark correlations at high densities providing a non- 
discontinuity in the physical observables, unlike the original version of the QMC model with 
constant bag. This suggests that in a QMC model with bag constant coupled to the scalar 
field, one can in principle extrapolate to high densities to explore the neutron star properties. 



IV. COMPOSITION AND STRUCTURE OF NEUTRON STAR MATTER 

In this section the constitution and structure of stable charge neutral (neutron star) 
matter in the supernuclear density regime are presented in the QMC and QHD models. In 
Fig. 6 the abundances of baryons and leptons as a function of density in the star matter are 
shown for the models QMCII (top panel) and QHDII (bottom panel). At densities slightly 
below the nuclear matter value, the /?-decay of neutrons to muons are allowed, and thus 
muons start to populate. The charge neutrality of a star forces a high isospin asymmetry so 
both the electron chemical potential /i e and the p field, p 3 ; grow at low density as evident 
from Fig. 5. Although the p field is very small and never exceeds —18 MeV, its correct 
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determination from symmetry energy is of utmost importance as it determines the proton 



fraction. It was demonstrated |41|] that for a npe system rapid cooling by nucleon direct 
URCA process is allowed by the momentum conservation condition k p + k e > k n which 
corresponds to a proton fraction Y p > 0.11. In both the models used here this condition is 
satisfied at densities ub > 0.28 fm~ 3 thus rapid cooling by direct URCA process can occur. 
In the absence of any hyperons, the charge neutrality condition forces the proton (and lepton) 
fractions to continuously increase with density in these relativistic mean field models. Thus 
once the threshold density for cooling by direct URCA process is achieved, it would persist 
up to the center of such stars. In contrast, the decreasing symmetry energy at high densities 
in the nonrelativistic models would limit cooling near the central region of massive stars. The 
symmetry energy also determines the hyperon production threshold density obtained by the 
condition /i# = fi n — <lB^e > Ssik = 0), where the energy eigenvalue Eb is given by Eq. (24). 
Consequently the threshold density for a hyperon species is determined by its charge and 
effective mass and by all the fields present in the system. As expected, the A with mass 1116 
MeV and S~ with a mass 1193 MeV appear at roughly the same density % ~ 0.38 fm -3 , 
because the somewhat larger mass of £~ is compensated by its negative charge. Since charge 
neutrality can now occur more economically by £~, the lepton fraction begins to fall. The 
electron chemical potential (see Fig. 5) then saturates around 200 MeV and subsequently 
decreases with increasing £~ population. More massive and positively charged particles than 
these appear at high densities. The substantial reduction in the effective mass of S in QHDII 
(as shown in Fig. 4) is manifested by a relatively early appearance of H~ and S° and their 
larger abundances in the star compared to that in QMCII. The enhanced S~ production 
in turn causes a further rapid decrease of the lepton fraction in QHDII. At high densities, 
all the baryons tend to saturate with the abundance of A being maximum, even exceeding 
the number of neutrons. Because of the fast growth of the hyperons and their comparable 
abundances with the nucleons, the dense interior of a star resembles more a hyperon star 
than a neutron star. The net strangeness fraction for stars in the model QHDII is slightly 
enhanced due to larger H abundance than that in QMCII. In contrast to pure npe stars, the 
proton fraction here reaches maximum value once the hyperons (E, A) start populating and 
thereafter it saturates at the level of 20%. Therefore rapid cooling by nucleon direct URCA 
process can still occur in these stars. Since the critical density of nucleon direct URCA 
process is nearly identical to the hyperon threshold density, and the emmisivities from the 
hyperon direct URCA processes are about 5 — 100 times smaller than that from the nucleons 
pE2]| , the stars cooling by direct URCA process is dominated by nucleons — the hyperons 
would only have a minor contribution to it. In the models QMCI and QHDI, the effective 
masses of the hyperons and the potentials are however found to be almost identical, hence 
the composition of the stars in model I were found to be practically indistinguishable. 

The equation of state, pressure P versus the energy density e is displayed in Fig. 7 for 
the different models studied here. The EOS for nucleons only (np) star shown for the QMC 
model (solid line) is found to be considerably stiff. Since the corresponding EOS in the 
QHD model is found to be nearly identical we do not present the result for clarity. At high 
densities when the Fermi energy of nucleons exceeds the effective mass of hyperons minus 
their associated interaction energy, the conversion of nucleons to hyperons is energetically 
favorable. Since this conversion relieves the Fermi pressure of the nucleons, the equation of 
state is softened. This effect is further accentuated by the decrease of the pressure exerted by 
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leptons because of their replacement by negatively charged hyperons in maintaining charge 
neutrality more efficiently The EOS for nucleons plus hyperons (npH) system for both the 
models I and II are shown in the figure; the thick lines refer to QMC results while the thin 
lines correspond to QHD. Several structures observed in the equation of state correspond to 
the densities at which different hyperon species begins to populate. In model I (shown by 
dashed lines), the EOS in QMCI is found to be softer compared to QHDI. In this situation 
the vector fields (u, p) and the effective hyperon masses are found to be same in the two 
models, only the scalar field a and the effective mass of nucleons in QMCI is somewhat larger 
than QHDI. This results in the QMCI a larger contribution from the scalar attraction and 
a smaller one from the repulsive kinetic term of nucleons to the pressure (Eq. (26)) leading 
to a softened EOS. With the inclusion of strange mesons, the effective masses of hyperons 
also undergo significant reduction (see Fig. 1) while, in general, the contribution from the 
repulsive field dominates over the attraction from a* field. The net effect is thus a stiffer 
EOS in model II (shown by dash-dotted lines) in contrast to model I. As evident from Figs. 
4 and 5, the combined effects of enhanced effective masses, considerably large repulsive 
field and smaller attraction from a* field act in increasing the pressure in QMCII compared 
to QHDII. This entails a pronounced stiffening with the EOS for QMCII being even stiffer 
relative to QHDII. Thus by including the strange mesons, a complete reversal in behavior 
occurs for the equation of state in QMC and QHD models which should have a significant 
bearing on the structure of the stars. For comparison, the causal limit p = e is also shown 
in the figure. All the relativistic models studies here respect causality condition dp/de < 1 
so that the speed of sound remains lower than the speed of light. 

An important parameter describing the equations of state is the adiabatic index T = 
dlnP/dlmiB = (P + e)/P ■ dP/de. In Newtonian theory it is possible to find a stable 
hydrostatic configuration for a spherical mass distribution if the adiabatic index exceeds 
4/3 |43| , and when general relativity is included it slightly exceeds this value [44]. Fig. 8 



shows the adiabatic index Y versus the energy density e for different equations of state. At 
densities below the nuclear matter value, T could have very small values since most the 
pressure support for the star originates from the electrons, the nuclear phase has a negative 
contribution to pressure. At densities greater than about n , the EOS stiffens and T is 
significantly greater than 4/3. With the appearance of hyperons, the softening of the EOS 
is manifested by the considerable lowering of the adiabatic index. The several structures 
observed in the EOS corresponding to the population of hyperons at the threshold densities 
is clearly evident in Fig. 8. The adiabatic index drops at each density when a new hyperon 
species is populated. 

The differences in the EOS at high densities are expected to be reflected in the structure 
of the neutron stars, namely their masses and radii. The static neutron star sequences 
representing the stellar masses M/M & and the corresponding central energy densities e 



c 



obtained by solving the Tolman-Oppenheimer-Volkoff equations |45| are shown in Fig. 9 
for different equations of state. In general, such a sequence posses a minimum mass below 
which gravitational attraction is not sufficient against the radial oscillations that destroy 
these configurations by dispersal. On the other hand, a maximum mass of the sequence 
exits beyond which the pressure support from the EOS is insufficient against the strong 
gravitational attraction. Stars beyond this mass are unstable to acoustical radial vibrations 
and thereby collapse to a black hole. The crustal region has a negligible contribution (~ 
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1(T 5 M ) to the total mass of a star, while most of the mass originates from the dense 
interior beyond the saturation density. Thus mass measurements may provide considerable 
insight into the interior constitution of a star. For the np system, the extremely stiff EOS 
corresponds to a large Fermi pressure and hence can sustain large limiting mass. The 
maximum masses for such np stars in the QMC model is M max = 1.988M , while a relatively 
softer EOS in the corresponding QHD model results in M max = 1.962M©. In Table III, the 
maximum masses and the corresponding radii i?jw max and central baryonic densities n c are 
presented. In Fig. 10, we also show the mass-radius relationship of the different EOS. With 
the inclusion of more baryon species in the form of hyperons, the considerable softening of the 
EOS results in relatively much smaller mass stars. Since the QMCI model has a much softer 
EOS than QHDI, the M max values are 1.478M© and 1.488M respectively; the star sequence 
and the mass-radius relationship in these models are not shown in the figures for clarity. The 
larger radii obtained in stars with hyperonization are a consequence of weaker gravitational 
attraction from the smaller masses that causes the stars to be large and diffuse. Compared 
to the mass of a star where the contribution is primarily from beyond the saturation density, 
about 40% of the star's radius originates from the EOS at ub ~ n . Consequently, radius 
measurements should be rather insensitive to the changes in the EOS due to quark structure 
of baryons or to the interior constitution of the star. In fact, no precise radius measurements 
currently exist. The star sequences and the mass-radius relationships in model II are shown 
in Figs. 9 and 10, respectively. The reversal in behavior observed in the EOS with the 
addition of strange mesons, i.e. the EOS is stiffer in QMCII than QHDII, is manifested in 
the maximum masses and the corresponding radii for such stars. For all the cases studied 
here, the maximum masses of the stars are found to be larger than the current observational 



lower limit of 1.44M imposed by the larger mass of the binary pulsar PSR 1913 + 16 |46 



When hyperons are included, the central densities reached for all the stars above this lower 
mass limit suggest that the inner cores of these stars are rich in hyperons. 

Constraints on the high density EOS can be derived from the measurement of the absolute 
upper limit to the rotation velocity f2 max of a neutron star. This is possible because the 
maximum angular velocity of rotating neutron stars, in general, is an increasing function of 
the softness of the EOS which corresponds to more compact stars with high central densities 
and smaller equatorial radius. The maximum rotation rate of a neutron star is determined by 
the condition that the equatorial surface velocity equals the Keplerian velocity — the orbital 
velocity of a particle at the equator. At the Keplerian frequency, the rotating star is unstable 
with respect to mass shedding from the equator. In reality, the rotation may be more 



severely limited by gravitational radiation instability to nonaxisymmetric perturbations [47 



However, this instability has been later shown to be stabilized by the existence of the viscosity 
of stars at homogeneous density fl8"|| . Therefore, the Keplerian rate may be considered as 
a reasonable estimate of the maximum rotation rate of a neutron star. The calculation 
of the Keplerian velocity for a given EOS is quite involved and has to be performed with 
full general relativity. On the other hand, a precise universal empirical formula was found 
by Haensel and Zdunik [45] for the maximum angular frequency f2 max for rigid rotation 
in terms of maximum mass and radius of a nonrotating star for a given EOS: f2 max = 
Cn(M max /M & y^ 2 (RM m!iX /lO km)~ 3 / 2 s _1 , where M max and -RA/ max are the maximum mass and 
corresponding radius of the nonrotating star. The dimensionless phenomenological constant 
independent of the EOS was found P| to be Cq = 7750. Later calculations within the full 
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framework of general relativity have shown [50| that this formula has an accuracy better 



than 5% for a wide range of realistic EOS which are both causal and stiff enough to support 
observed maximum mass M max = 1.44M Q . Since rotational energy stabilize a star, the 
most massive rotating star has more mass and radius than the maximum mass and radius 
of the corresponding nonrotating star. For our present equations of state, the Keplerian 
frequencies obtained by using the above formula are given in Table III for the respective 
rotating stars. It is observed that f2 max closely follow the trend of the masses (and radii) of 
different models with the nucleons only star having the largest rotation rate and stars with 
hyperons in model I are the slowest. As mentioned above, this may be attributed to the 
smallest mass and thereby largest radius due to small gravitational attraction in the softest 
equation of state and vice versa. 



V. SUMMARY AND CONCLUSION 



In this paper we have investigated the effects of the internal quark structure of baryons on 
the neutron star properties within the relativistic mean-field quark-meson coupling model. 
This model describes baryons as nonoverlapping MIT bags in which the quarks interact 
through scalar and vector mean fields. Before any reliable prediction of the quark structure 
effects on neutron star properties at high densities can be made, it is essential that the QMC 
model reproduces at around nuclear matter densities the results of more established quantum 
hadrodynamics model where the relevant degrees of freedom are the structureless baryons. In 
the original QMC model, the bag constant was fixed at the free space value as a consequence 
of which the model predicts much smaller attractive scalar potential and hence smaller spin- 
orbit potential compared to the experimental results and that obtained in the QHD model. 
By considering a medium (density) dependent bag constant parametrized through direct 
coupling to the scalar field, correct spin-orbit splitting was observed. Moreover, this medium 
modified QMC model is found to be in excellent agreement with the low and moderate 
density results of QHD model with a general nonlinear scalar potential when both the models 
are calibrated to produce the same nuclear matter saturation properties. This improved fit 
is obtained by employing a bag constant which is a decreasing function of density, which 
however implies an increasing bag radius. We are therefore faced with the problem that 
due to the increasing bag radius and thereby overlapping bags the QMC model may not be 
applicable at the high density regime relevant to central densities of massive stars. A natural 
test of the reliability of the model at high densities then lies in the fact that the equation of 
state should be well-behaved and continuous when extrapolated to the extremes of density. 
Indeed, we have found that the physical observables exhibit reasonable behavior without 
any discontinuity at high densities up to ps 10n studied here. The deviation of the 
results observed at high densities in the two effective field theoretical models, the QMC and 
QHD, with different underlying basic constituents may be interpreted as primarily arising 
from the crucial effects of the quark structure. The original version of the QMC model, 
where the important effect of quark-quark correlation associated with overlapping bags is 
neglected, exhibited discontinuity and therefore could not be extrapolated to high densities. 
A direct coupling of the bag constant to the scalar field is equivalent to higher powers of 
(nonlinear) scalar interactions. According to the modern viewpoint of effective field theory, 
higher order meson interaction includes the effect of correlations. Therefore, by employing 
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a medium (scalar field) dependent bag constant, we not only could reproduce the correct 
scalar potential but could also mimic the quark-quark correlations leading to a well-behaved 
EOS at high densities. 

We have included two additional (hidden) strange mesons which couple only to strange 
quarks in a baryon bag in the QMC model and to hyperons in the QHD model. The 
rather strong hyperon-hyperon interaction can be accounted by these mesons. The coupling 
constants of the quarks and hyperons have been fixed by SU(6) symmetry relations based 
on quark-counting argument. The strange mesons are found to have considerable influence 
on the composition and structure of neutron star matter with hyperons. In absence of these 
mesons, the model QMCI exhibits softer EOS with smaller maximum mass star and larger 
corresponding radii than that in QHDI. With the inclusion of strange mesons, the additional 
attraction imparted by the scalar meson o* causes a drop in the effective masses of hyperons, 
the decrease being determined by the number of strange quarks in the baryons. However, 
the two mesons (a*, 0) has an overall repulsive effect so that the EOS in the model II 
is stiffer compared to model I without the strange mesons. The repulsion is maximum in 
QMCII because of the decreasing in-medium scalar-baryon coupling constant and smaller 
scalar field a* than in QHDII. Consequently, the EOS in QMCII is significantly stiffened 
(and even stiffer than QHDI) with relatively larger maximum mass star and corresponding 
smaller radius. As observed in previous studies, the EOS is found to be considerably softened 
by incorporating hyperons as the new degrees of freedom which appear at ~ 2n . The 
center of massive stars are found with comparable abundance of hyperons and nucleons; 
the strangeness fraction of stars in QMC models are relatively higher than in QHD models 
because of enhanced S production. Rapid cooling by direct URCA process of all these 
stars are found to be dominated by nucleons due to large proton fraction (Y p > 0.15), the 
hyperons add only ~ 5% to this more dominant process. It therefore seems difficult to 
differentiate stars in the different models studied here with and without hyperons by rapid 
cooling procedure. 
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TABLES 



TABLE I. The free space values of bag parameters zq and bag radii Ro for different baryons 
obtained by reproducing the baryonic masses mg in free space. The bag constant for the baryons at 
the free space value is B^ 4 = 188.1 MeV while the masses of the quarks are taken as m u = mid = 



and m s 


= 150 MeV. The strangeness 


Sb of the baryons are also 


given. 




Bar yon 


m B (MeV) 


S B 


zo 


R (fm) 


N 


939 





2.030 


0.600 


A 


1116 


-1 


1.815 


0.642 


£ 


1193 


-1 


1.629 


0.669 




1313 


-2 


1.505 


0.686 



TABLE II. The coupling constants obtained in the QMC model by reproducing the saturation 
density uq = 0.17 fm -3 , the binding energy B/A = 16 MeV and the symmetry energy a sym = 33.2 
MeV. The coupling constant of the scalar a field to the bag is = 2.269. The scalar coupling 
constant corresponds to the free space value [see Eq. (20)] while its coupling to the («, d) quarks is 
taken as g% = 1. The predicted values of compressibility and effective nucleon mass at saturation 
are K = 289 MeV and m* N /mN = 0.78. The coupling constants in the QHD model are obtained 
by adjusting these same five saturation properties. The meson masses are taken to be m a = 550 
MeV, m w = 783 MeV and m p = 770 MeV. 



Model 


9l/^ 


SS/47T 




92 (fm l ) 


53 


QMC 


5.184 


5.240 


5.203 






QHD 


5.174 


5.339 


5.146 


12.139 


48.414 



TABLE III. The maximum masses M max /M0 of nonrotating stars and their correspond- 
ing radii i?M max and central densities n c in the models QMC and QHD. The Keplerian 
frequency fl max for the respective rotating configurations are obtained from the relation 
fi max = 7750(M max /M Q ) 1 / 2 (i?M max /10 km) -3 / 2 s" 1 . Results are for stars with only nucleons (np); 
stars with further inclusion of hyperons (npH) in the model I where the interaction is mediated by 
(cr, u, p) mesons, and stars in model II where the additional mesons (a*, (f>) are included. 





QMC 


QHD 




M max /M 




n c 




M max /M 




n c 


^max 






(km) 


(fm" 3 ) 


(10 3 s- 1 ) 




(km) 


(fm" 3 ) 


(10 3 s- 1 ) 


np 


1.988 


10.632 


1.102 


9.969 


1.962 


10.561 


1.110 


10.001 


npH (model I) 


1.478 


11.242 


0.965 


7.904 


1.488 


11.092 


0.999 


8.091 


npH (model II) 


1.539 


10.823 


1.096 


8.538 


1.491 


11.040 


1.022 


8.159 



22 




rn E /m s 

m^/m s 
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FIG. 1. The variation of effective nfa'sse^mgy m^ 6 of the barybhs a$ a function of baryon density 
in the models QMCI (thin lines) and QMCII ^Ihick lines). 




FIG. 2. The variation of bag constant Wb/Bq of"the°rjar^ons 1 fe a 1 function of baryon density ns 
in the models QMCI (thin lines) and QMCII (thicl! fines)'. The free space bag constant is taken as 



B, 



1/4 



188.1 MeV. 
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FIG. 3. The variation of bag radiu§"°i?£ /Ko°6 4 i' the 1 bSfvpris 1 as 1 a fhrfetion of baryon density in 
the models QMCI (thin lines) and QMCII (thick fines?. The free space bag radius Rq for different 
baryons are given in Table I. 
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0.5 
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FIG. 4. The variation of effective masse¥ m B 4 / m°B ofTlje rJarybhs as a function of baryon density 
us in the models QMCII (thick lines) and QHDII B (fhin lines). 
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FIG. 5. The mean meson fields and the electrochemical potential versus the baryon density tib 
for the models QMCII (top panel) and QHDII (bottom panel). 



25 




26 



np (QMC) 

npH (model 1) 

npH (model II) 
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FIG. 7. The equation of state, pres§ure P f§. energf density for Mfcleons only (np) star (solid 
line) in the QMC model, and for stars with further inclusion of hyperons (npH) in models I (dashed 
lines) and in models II (dash-dotted) lines. The results are for QMC models (thick lines) and QHD 
models (thin lines). The causal limit (P = e) is also shown. 




FIG. 8. The adiabatic index T vSj the energy density e.for nucleonsl only (np) star (solid lines), 



U OUU I UUU IOUU ilUUU 

and for stars with further inclusion of hyperons ^npH) in yodels I (dashed lines) and in models II 
(dash-dotted) lines. The results are for QMC models (thick lines) and QHD models (thin lines). 
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FIG. 9. The neutron star sequence's^ tie 
and for stars with further inclusions,^ 




npH (model II) 



the limiting mass for nucleor.s only (np) star (solid lines) 
igpOKgng^npH) ^th& n&dejJ i>[ (dash-dotted lines). The 



results are for the models QMC (thick lines) aftd'Qffl^tliin lines). The filled circles correspond 
to the maximum masses. 



FIG. 10. The mass-radius relation for neutron star sequences near th 




npH (model II) 



e limiting mass for nucleons 



only (np) star (solid lines) and for°s, tars with further, inclusion of hyg grons (npH) in the model II 
(dash-dotted lines). The results are for the models T^KlC (thick lines) and QHD (thin lines). The 
filled circles represent the maximum masses. 
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